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We explore the quantum phase transitions between two ordered states in the infinite dimensional Hubbard-Holstein 
model at half filling. Our study is based on the dynamical mean field theory (DMFT) combined with the numerical 
renormalization group (NRG), which allows us to handle both strong electron-electron and strong electron-phonon 
interactions. The transition line is characterized by an effective electron-electron interaction. Depending on this effec- 
tive interaction and the phonon frequency uiq one finds either a continuous transition or discontinuous transition. Here, 
the analysis focuses on the behavior of the system when the electron-electron repulsion U and the phonon-mediated 
attraction A are equal. We first discuss the adiabatic and antiadiabatic limiting cases. For finite luq we study the differ- 
ences between the antiferromagnetic (AFM) and charge order, and find that when present the AFM state has a lower 
energy on the line. 
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1 Introduction A feature of strongly correlated sys- 
tems is the existence of competing interactions on low en- 
ergy scales which can lead to different types of symmetry 
breaking and different ground states. There can be transi- 
tions to various forms of magnetic order, to superconduct- 
ing or charge ordered states; there also may be transitions 
between these states, and in some cases they even coexist. 
Here we study the competing effects arising from the inter- 
electron interactions and the electronic coupling to lattice 
modes, as described by the Hubbard-Holstein (HH) model, 
and the competition between antiferromagnetic (AFM) and 
charge order (CO) at half filling. Our study is based on the 
dynamical mean field theory (DMFT) combined with the 
numerical renormalization group fl]|2] (NRG). The DMFT 
becomes exact in the limit of infinite dimensions [3|, and it 
can generate non-perturbative solutions, such that electron- 
electron and electron-phonon interactions with arbitrary 
coupling strengths can be studied. 

There have been several applications of the DMFT- 
NRG method to study phase transitions in the Hubbard- 
Holstein model JUEEEM] ■ There are various possible 



transitions to states of broken symmetry in this model; 
bipolaronic (BP), CO, AFM and the superconducting (SC) 
state. We restrict our attention here to the case of half- 
filling and zero temperature. The transition first studied by 
the DMFT-NRG method for this model did not include the 
possibility of either CO or AFM |9,6|. There is, however, 
a metal-insulator transition from the normal (N) to the BP 
state [6|. If the possibility of transitions to CO and AFM 
are included, then it is found that as the attractive term in- 
duced by A = 2g 2 /uq overcomes the repulsion due to U 
the system changes the ground state from AFM to CO iflOl 
[TTI . Close inspection shows that the transition between the 
CO and the AFM state generally does not exactly occur 
at A = U, but at slightly larger values of A. Our results 
show evidence for a direct transition from an ordered to 
an ordered state. There can be a continuous transition for 
smaller values of the interactions U and g and discontin- 
uous transitions for larger interactions. The focus of this 
paper will be the behavior along the line U = X. 
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2 Model and method The Hamiltonian for the 
Hubbard-Holstein model is given by 



JT =-*Z)(cJ, ff C;> + ^cO 



U 
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c\ creates an electron at lattice site i, and b\ a phonon 

with oscillator frequency wo, n i.<r = cj a Ci^ a . The elec- 
trons interact locally with strength U, and their density 
is coupled to an optical phonon mode with coupling 
constant g. We assume a semi-elliptic DOS, po(s) = 
2y/D 2 - e 2 /nD 2 , with D = 2t. t = 1 sets the energy 
scale in the following. 

For our calculations we assume a bipartite lattice with 
A and B sublattice, where the matrix Green's function can 
be written in the form 



GkJoj): 
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(A,<r(u)(B,a{u) 






(2) 

With ( a ,cr(u) = U) + /ia,er — -S'o:,cr(w), OC = A, B, and k- 

independent self-energy. For commensurate charge order 
we have \ia,o = At — /i c , /i_B, CT = H + h c and .17B :Cr (u;) = 
t/n-.EU.o-C-w)*, with n = (n A +n B )/2, n a = J2 a n a,a- 
For the AFM order one has [iA,a = A* ~~ ^sj Ms.cr = 
/i+(t/i s , and the condition Z'b.ct (a;) = Sa.-u (^)- We con- 
sider solutions where the symmetry breaking fields vanish, 
h c , h s — > 0. In the AFM case the A-sublattice magnetiza- 
tion, <P a f m = rriA — {nA.\ — tia,i)/2 serves as an order 
parameter. For CO we define <P co = (iia — l)/2. 

In the DMFT this local Green's function, and the self- 
energy are identified with the corresponding quantities for 
an effective impurity model 0. One focuses for the cal- 
culations on the properties of the A-sublattice. We can 
take the form of this impurity model to correspond to an 
Anderson-Holstein impurity model |[T2l and calculations 
are carried out as detailed, for instance in Ref. [13, 14 1 . We 
solve the effective impurity problem with NRG adapted to 
these cases with symmetry breaking. For the logarithmic 
discretization parameter we take the value A = 1.8 and 
keep about 1000 states at each iteration. The initial bosonic 
Hilbert space is restricted to a maximum of 50 states. 

3 Behavior along U — A As the electron-phonon 
coupling in (Q]i is linear, the bosonic field can be integrated 
out in a path integral framework, which yields an effective 
electron-electron interaction of the form 



U eS {uj) = U + 



2g 2 tu 



(3) 



The two terms are the competing interactions on different 
energy scales. For large u>, the electron repulsion U dom- 
inates. Near u = u?q the retarded effective attraction due 



to the phonon starts to play a role and in the limit w — >• 
the expression tends to f/ e ff = U — A. We are studying 
the case U — A, so U e g — at low energy. For the gen- 
eral behavior on all energy scales, finite uiq, and arbitrary 
coupling strength no simple description is available. First 
insight can be gained by considering limiting cases. 

In the antiadiabatic limit, luq —> 00, where A = 2g 2 /wo 
is kept fixed, U e g(w) becomes completely independent of 
w and tends to U c s = U — A, so that the model then be- 
comes equivalent to a pure Hubbard model with U = U e g. 
The system is then not ordered for any value U — A. For 
any finite U c s > the system is AFM ordered ifTBI and for 
Ueg < in the CO or SC state [[6). We have therefore a 
continuous transition from an ordered to an ordered state. 
In the DMFT-NRG calculations we find that the smaller 
wo is the larger the order parameter near the transition be- 
comes. Thus we conclude that this transition scenario per- 
sists for weak coupling and finite uiq. 

Another limit the adiabatic limit, luq —> 0, can be stud- 
ied transparently in static mean field theory. The mean field 
self-energy reads, 



I7 a ,,(w) = Un% - \{n a - 1), 



(4) 



independent of w. In order to determine n™ = (n"), we use 
the self-consistency equation 
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We have used the spectral function of the first element of 
the matrix Green's function for the bipartite lattice in the 
form (|2), 



Ca,a(^ + )/7T 



Pa,kAv) = - In V,.r,.,+v„_r,.rn- g 2 . < 6 > 



with lu + = lu + ir\. We have introduced 

£A,cr — ^A,cr + £b,<j — HB,c 



where 



±E a {e k ), 
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The expression for the total energy reads at half filling 



(8) 
(9) 



EnJ = E£l + (U- 2\)<P 2 co - U$l (m + ^, 
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where E^ is the kinetic energy term fill - We study ex- 
clusive order here, i.e., either CO or AFM; coexistence is 
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energetically not favorable. We can see that for order pa- 
rameters of equal magnitude the CO state has lower en- 
ergy for X > U and the AFM state otherwise. For U = X 
both states are degenerate. When we increase the interac- 
tion U = X from zero we find that the order increases as 
illustrated in Fig.Q] 




Figure 1 (Color online) The mean field expectation values 
•P along the U = A-line. Inset: Total energy. 

We see the typical mean-field exponential increase. One 
can, however, note a suppression due to the competing in- 
teraction present. The total energy is shown as an inset. 
The static mean field solution in the adiabatic limit gives 
always a discontinuous transition from CO to AFM state 
when U or A are changed. 

We now turn to the situation for finite ujq. This situa- 
tion is accessed with DMFT-NRG. In Fig. [2] we show the 
result for ujq = 0.2. For small values of U = A both order 
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Figure 2 (Color online) The expectation values <P along 
the U = A-line for ujq = 0.2. Inset: Total energy. 

parameters are zero. From a certain value of the interac- 



tion strength on solutions with finite AFM order can be 
found, and for even larger interactions also solutions with 
finite CO exist. The point in the phase diagram separat- 
ing first and second order transition lies close to this point 
where both order parameters are finite, Af s ~ 2.7. It will be 
shifted to larger values Af s when uq increases. Our numer- 
ical results give a continuous rise for ^ a f m , and a sudden 
increase for <£ co . 

From the total energies (inset) we can see that the AFM 
state has lower energy and is therefore favorable along the 
U = A-line. Calculations in the normal state reveal renor- 
malized quasiparticles (z < 1) with a small repulsive ef- 
fective quasiparticle interaction, which is consistent with 
this observation ifTTIl . Thus, there is an asymmetry of AFM 
and CO state for finite uq . 

In summary we have analyzed the competing interac- 
tions in the HH model and focused on the situation where 
U = X. For u?q — > oo there is no order along this line in 
the U — A-phase diagram, whereas for wj^O and finite 
interactions the system is ordered with either AFM or CO, 
which are degenerate. For intermediate values of luq there 
is no order for small couplings and AFM or CO can be 
present for larger couplings, where the AFM state is found 
to have lower energy. 
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